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<N Abstract 

A two-dimensional cellular automaton is introduced to model the flow and 
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q jamming of vehicular traffic in cities. Each site of the automaton represents 
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ON a crossing where a finite number of cars can wait approaching the crossing 

"bib 

from each of the four directions. The flow of cars obeys realistic traffic rules. 
Q_i We investigate the dependence of the average velocity of cars on the global 

a 

traffic density. At a critical threshold for the density the average velocity 

reduces drastically caused by jamming. For the low density regime we provide 

analytical results which agree with the numerical results. 
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I. INTRODUCTION 



Vehicular traffic on highways as well as in cities tends to suffer from a jamming transition 
when the global traffic density exceeds a critical threshold value. The phenomena related to 
traffic jams have attracted the attention of engineers and physicists since many years and 
there exists a large variety of experimental observations, see [1-4] and references therein. 
The formation of traffic jams has been investigated by many authors using various methods. 
As early as 1955 Lighthill and Whitham [5,6] described the spontaneous formation of regions 
of increased car concentrations and shock waves using their approach of kinematic waves, 
i.e. wave motion where the spatial value of the flow is a function of the spatial concent- 
ration distribution. Prigogine and Herman [7,8] investigated traffic jams on highways by 
performing a stability analysis of hydrodynamic equations. One-dimensional cellular au- 
tomata models [9] for the simulation of traffic flow have been proposed e.g. in [10-12]. 
This type of models has been studied intensively using massively parallel computers [13-17]. 
Schadschneider and Schreckenberg [18] solved a one-dimensional cellular automaton mo- 
del analytically in the mean field approximation, Csahok and Vicsek [19] investigated their 
automaton in the presence of quenched noise. Nonlinear wave descriptions have been pro- 
posed in [20-22]. Ben-Nairn et. al [23] apply a ballistic aggregation process to model the 
kinetics of clustering in one-dimensional traffic flow. Hydrodynamic approaches have been 
investigated by various authors, e.g. [24,25]. These models are based on certain assumptions 
concerning the "hydrodynamic" properties of traffic flows such as the velocity versus density 
relation [26,27], and the viscous terms [22]. 

One-dimensional traffic has been studied extensively and the developed models are suffi- 
ciently sophisticated to reproduce experimental observations, i.e. the fundamental diagram 
(throughput versus flow density) of one lane roads. Some of the authors (e.g. [28]) claim 
that the size distribution of the jams obeys a power law, but very large scale computations 
have shown that the lifetimes of the jams and hence their sizes reveal a characteristic cut- 
off which seems to be no finite size effect of the simulation [14]. There are several models 
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describing traffic flow in more than one dimension. The analysis of traffic on a network with 
multiple sources and sinks can be found e.g. in [29-32]. Nagatani [33] describes a cellular 
automaton representing a two-lane roadway. 

An interesting new approach was currently proposed by Bando et al. [34] who express the 
rules which govern the behavior of the cars (i.e. acceleration and deceleration) by a dynamic 
equation for each car x n = a {V (x n+ \ — x n ) — x n }. The interaction of neighboring cars is 
expressed in terms of the function V of the distance between the cars. Using this approach 
the authors connect in some respect the model of the cellular automata with molecular 
dynamics. 

Two-dimensional cellular automata designed to simulate traffic in a city have been inves- 
tigated by some authors. Biham et al. [35] proposed a three state cellular automaton model 
where at a given time step each site can be occupied by a car moving from South to North, 
by a car moving from West to East, or the site may be empty. The time behavior is ruled by 
synchronous traffic lights at each site, allowing alternatively for vertical or horizontal traffic. 
Obviously there are situations when the traffic jams, i.e. when the traffic light allows for 
driving but the next site is occupied. The model given in [35] and particularly the spatial 
correlations in the jamming phase have been numerically investigated by Tadaki and Kiku- 
chi [36]. A similar model but with "faulty" traffic lights was studied in [37]. Nagatani [38] 
investigated the spreading of a jam which is induced by an accident using an extremely 
simple automaton rule. Numerically he finds scaling laws for the size of the spreading jam 
as a function of time elapsed since the occurrence of the accident which cannot be derived 
analytically so far. Fukui et al. [39] numerically investigated the evolution of ensemble ave- 
rages of the jamming process for the simple automaton model described in [35]. For the 
same model [35] Chau et al. [40] recently gave an analytic upper bound (depending on the 
dimension) for the critical car density r] cr , i.e. when the system transits from the "moving 
phase" into the "jamming phase". In two dimensions they found r\ cr < 11/12. Obviously 
this upper bound is not close to the value which was observed in numerical simulations. In 
the case of Nagatani and Seno [41] there is a set of parallel one-way streets, all oriented 
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in the same direction (rr-direction), which is intersected by a single perpendicular one-way 
street in ^-direction. 

For fixed car density on the perpendicular street p y they find that the flux of cars J x 
rises linearly with the car density p x until a characteristic critical threshold p x is reached. 
When further increasing the density p x the flux J x approximately remains constant while 
the average velocity of the cars (v x ) drops. This rather sharp change in the overall behavior 
is due to the formation of traffic jams where the flow is irregular and discontinuous. When 
further increasing the density p x the discontinuous character of the flow disappears at a 
second threshold p^ and the flux J x declines linearly. This behavior corresponds to the 
symmetry of the car density in ^-direction and the spaces between the cars 1 — p y . Indeed, 
one finds that the critical densities p c x and p x are almost exactly symmetric, e.g. for p y = 0.3 
one finds p x = 0.31 and p x = 0.69. 

Cuesta et al. [42,43] have been the first who reported on simulations of a two-dimensional 
automaton where the cars are allowed to change their direction. In their model each crossing 
can be occupied by one car or it can be empty. Each car is assigned a preferred direction; 
the parameter Wi(r) is the probability that the i-th car at the position r moves in the 
next time step horizontally and 1 — Wi(r) is the probability to move in vertical direction. 
Horizontal motion is allowed at even time steps, vertical motion at odd times. There have 
been defined two variants, A: there are only one-way streets directed from South to North 
and from East to West. Half of the cars is given the trend Wi(r) = 7 and the other half is 
given Wi(r) = 1 — 7. In the second variant B there are one-way streets pointing alternatively 
South — > North and North — > South, and West — > East and East — > West, respectively. The 
trends of the cars which are subdivided into four equal groups, point into one of the four 
possible directions. 

We should mention that several authors claim that at least the one-dimensional traffic 
flow problem is closely related to sand flowing in pipes, e.g. [44-48]. 

In the present paper we are concerned with a two-dimensional cellular automaton desi- 
gned to simulate traffic in a city. As a special feature of our model we emphasize the fact 
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that we use rather realistic traffic rules which will cause a slowing down of the average traffic 
velocity and finally will lead to the collapse. 

II. THE MODEL 

In our model the city is represented by a set of L streets in horizontal direction crossing 
L streets in vertical direction. Cars are allowed to move in both directions, i.e. we have no 
restrictions to one-way streets. The crossings define a two-dimensional cellular automaton, 
each of the L 2 crossings is represented by a site of the automaton. We assume periodic 
boundary conditions in both directions. Fig. 1 shows a schematic plot of an automaton site 
(a crossing). At each crossing there are four queues of maximum length Q filled by cars 
coming from one of the four directions, representing the finite space where cars can move 
freely between crossings in realistic urban traffic systems. Hence, at each crossing there are 
allowed at most 4 Q cars. Each of the cars has a desired direction, i.e. to the right, to the 
left or straight on, which has to be chosen according to a certain rule. We will discuss this 
point in detail below. Within each time step the first of the cars in each queue can go to the 
next crossing provided that it does not have to give way to one of the other cars at the same 
site. In case it has to give way to another car it will not move during this time step. In our 
model we assume the simple and realistic rule that a car has to give way in the case that 
another car occupies the same crossing at its right hand side. Moreover a car has to stop 
if it intends to turn left and there is a car at the opposite side of the crossing which goes 
straight on or turns to the right. If all of the four top positions of a crossing are occupied 
by four cars one of the cars will be selected randomly and will then be allowed to move in 
the current time step. Insofar we have chosen realistic rules which are current law in many 
countries. 

Each of the N cars in the system starts at a randomly selected crossing, its desired 
direction (left, straight on, right) will be determined according to some rule (discussed 
below). Within each time step a car maximally moves a distance of one lattice cell, i.e. it 
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tries to transit to one of the next neighbors of the current site provided the following three 
conditions hold: 

1. It does not have to give way. 

2. It is placed at the top position of its queue. 

3. There are less than Q cars standing in the destination queue. 

Otherwise it rests. Obviously, if there is only one car in town it never has to give way nor 
will it find another car in front of itself. Therefore the average velocity 1 (v) = 1 in units of 
lattice cells per time steps. 

As mentioned above we have to discuss the rule which determines the desired direction 
within each time step. There are at least three reasonable rules: 

1. The desired direction for each car is selected at random within each time step not 
regarding whether the car moved in the previous time step or stopped. This rule is 
rather artificial since provided a jam occurs it can dissolve due to the fact that the 
cars might give up their previously desired direction when they cannot move. 

2. Each car is assigned a destination site on the lattice. Once a car reached its destination 
it will be assigned a new randomly chosen destination. This rule is problematical since 
one has to define the detailed path a car follows to reach its destination. A natural 
choice would be to chose a path so that the car has to change its direction only once 
on its way between its starting position and its destination. In this case, however, 
the overall behavior of the system would depend on the size of the lattice L since one 
can check that the probability to be allowed to move straight on is higher than to be 
allowed to turn. Cars which intend to turn either left or right (on the average) have 



1 Here and in the following we describe the ensemble average of a value x by the symbol (x), 
whereas the time average by x. 
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to give way more frequently than those which go straight on. The larger the system 
the smaller becomes the probability per time step to turn. For the case L — > oo the 
probability for a car in a given time step to turn approaches zero. Hence this rule does 
not allow for size-independent results, and the results of a simulation using this rule 
apply only for one fixed lattice size L. 

3. Each car is assigned a desired direction during the initialization. But now, only after 
moving one step to the aimed destination a new desired direction will be chosen at 
random, i.e. the primarily chosen desired direction is maintained by each car until all of 
the three conditions hold which allow the car to move. This rule has several advantages 
compared with the other rules. First we avoid the size dependent behavior of the 
previous rule, and second we do not have long range correlations between the lattice 
sites for the case of low car density, which would be caused by the second rule. Hence 
the behavior of the system can be described by a Markov process, which will be the 
basic assumption of our analytical calculation in section IV. An analytical description 
which takes into account the long range correlation between the sites coming from the 
rules which determine the path of the cars seems to be very difficult, except perhaps 
for the (trivial) case L — > oo, when the probability of the cars to turn, and hence the 
traffic rules connected with the turn, can be neglected. In the following we always 
refer to the third rule. 

III. NUMERICAL RESULTS 

In this section we present the results found by numerical simulations of a cellular auto- 
maton which represents the traffic rules described above. 

The two-dimensional automaton consists of L x L sites on a rectangular lattice. To 
eliminate boundary induced effects we have chosen periodic boundary conditions in both 
directions, i.e. we assume the topology of a torus. We simulated systems of size L = 20, 
30, 50, 100 and no striking difference in the behavior of the system was observed. Hence, 

7 



we exclude finite size effects. The fact that already comparatively small lattices exhibit this 
independence from their size seems to indicate the existence of only weak spatial correlations 
between sites separated by a large distance. For high densities, i.e. close to or above the 
critical density r] cr , this statement does not remain valid. In fact, we expect the existence 
of long range correlations. However, in the present paper we are mostly concerned with 
the low density region and we will not devote ourselves extensively to the behavior of the 
automaton for high car density. 

Each of the L x L automaton sites had its own four queues of length Q = 10 connecting 
the site with one of its four neighboring sites as described in the previous section. For 
each value of the density 77 = j^j;, where N is the number of cars, we started the system 
at randomly chosen initial conditions. Each car was assigned the initial site, the queue in 
this site (i.e. whether it comes from North, East, South or West) and the desired direction 
(left, straight, right) for the next step. When a car transits to an empty site from a certain 
direction it occupies the first empty place in the queue. If it is at the top position of the 
queue it can transit to the next site or it has to wait, depending on the cars coming from 
the other directions of the same site, and depending on the traffic rules which apply to the 
given situation. If a car transits to another site and if there are other cars located in the 
previously mentioned queue, all of them advance by one position in this queue hence, moving 
up "closer to the crossing" . 

After starting the simulation for a given number of cars, i.e. for a given density 77, we 
evolved the system 10,000 time steps without recording the statistics in order to let the 
system forget about the initial preparation. Typically a "stationary" behavior was found 
after only a few hundred time steps. In density regions very close to the transition value 
however, the relaxation of the system lasts longer which hints at some phase-transition-like 
phenomenon. The mean velocity of the cars was then computed using 2,000 time steps. We 
varied the car density beginning from 77 = to 77 = 3.8 in steps A?7 = 0.01. 

The results of the simulations are shown in fig. 2. Like in simulations of cellular automata 
performed by other authors with different traffic rules we found a smoothly decaying function 
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in the low density region. At a certain threshold value rj cr ~ 1.6 the character of the traffic 
flow changes abruptly and the system transits into the clogged regime. But different from 
other systems we find a remnant velocity in the clogged regime. This is due to the rule that 
in case all four of the queues are occupied there will be one car selected randomly which then 
will be allowed to move. Therefore, in our model there exist situations where a jamming at 
a crossing with four queues occupied may still dissolve. 

The curve in the small density region is very smooth, i.e. there are almost no fluctuations. 
Hence, we can assume that long range correlations do not play a major role there. In the 
following section IV we will present a statistical description appropriate to model the flow 
in the low density regime, i.e. for i] <^ i] cr , which was inspired by this observation. 

IV. A STATISTICAL DESCRIPTION OF VEHICULAR TRAFFIC 

Let the number of cars be denoted by N and the number of cells respectively crossings 
building the torus by M = L x L where L is the common perimeter of the torus. For the 
global density we choose the symbol 77 = N/M. 

The quantity we focus our attention on is a time averaged 2 and normalized velocity 

1 f N-A(t) (A) T 

where A(t) means the total number of cars prevented from moving by traffic rules at time 
t which depends on the exact situation experienced during the simulation. 

For a description of the dynamics in the low density region r\ <C 1 we apply the following 
arguments: 

1. An effective description is possible which separates the process of clustering at the M 
possible crossings from the process of obeying the traffic rules. Due to this indepen- 



2 When performing the simulation the parameter T is chosen sufficiently large in order to scan the 
state space respectively the limit set with sufficient accuracy. 
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dence assumption the time average (1) extracted from a simulation can be decomposed 
into two averages: first a time and cluster average related to the traffic rules and second 
a time average with respect to the clustering process. 

2. Given i cars meeting at one crossing a variety of situations is possible. Every car is 
approaching the crossing from either North, East, South or West and intends to go 
straight, to turn left or to turn right. In the long run all possible situations will have 
equal probability. Hence, the average number of cars which must wait, denoted by 
Si with < S{ < i, can be computed by combinatorial reasoning. The result of this 
computation is collected in table II for i = 1, . . . , 8. A sketch of the derivation is given 
in appendix A. 

For clusters of size i > 8 the <5j are not shown since their contributions are suppressed 
by very small cluster probabilities. 

3. Central to our description is the identification of cluster distributions, i.e. we write 
k = (ho, hi, . . . , k N ) meaning that hi denotes the number of cells in the system where 
i cars are meeting. Of course, there are two boundary conditions, namely 

N 

E^ = M ( 2 ) 

i=0 

and 

N N 

J2ih = J2ik t = N . (3) 

i=0 i=l 

We now assume that for the low density regime i] <C 1 the process of clustering is 
ergodic which means for almost all initial configurations the time average can equiva- 
lently be represented by an average according to an invariant probability distribution. 
Note that this probability distribution is defined over the space of all configurations 
k compatible with the boundary conditions (2) and (3). In order to derive this distri- 
bution we additionally assume that the dynamics effectively is equivalent to a process 
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of independently distributing N cars (balls) among M cells with crossings (urns) ac- 
cording to the equidistribution. This means we assume the discrete dynamics to be a 
Bernoulli process. Then the probability distribution reads 

Ml Nl 

p(k) = ^ . (4) 

m n n m ki M 

To understand formula (4) in detail one has to realize that there exist 
Ml/ [k lkil . . . k N l] different ways to index the cells with numbers 0,1,..., iV since 
all cells containing the same number i of cars are not distinguishable. Then there 
exist Nl/ [(0!) fc °(l!) fcl . . . (AH) fcjv ] ways to index the cars with numbers 1, . . . , N taking 
care of the fact that cars within an identical cell are not distinguishable. Finally the 
equidistribution with respect to all configurations yields the factor M~ N . This result 
was obtained by von Mises (1939) [49]. 

Now the time average (A) T is replaced by the ensemble average 

* N 

{A) K = 5>(A;) 5>£«i- (5) 

k i=l 

The sum with respect to i reflects the average delay due to the traffic rules at the crossings 
and the sum with respect to all configurations obeying (2) and (3) - which is indicated 
by the symbol * - accounts for the cluster statistics. The factor accounts for the fact 
that the dynamics of the real process differs from a Bernoulli-process. In general there may 
exist spatial as well as temporal correlations. The spatial correlations become important 
for densely filled lattices. The dynamics which control the traffic within one cell is strongly 
influenced by the cars moving in neighboring cells. Finally the emergence of a traffic jam 
clearly expresses these correlations. On the contrary temporal correlations play the domi- 
nant role for very sparse system, i.e. in the low density region. Consecutive situations are 
not completely statistically independent. Let us briefly explain this effect considering only 
cluster of size i = 2: When the cars are randomly distributed amongst the lattice sites there 
is a certain probability that two cars p and q occupy the same site and hence form a cluster 
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of size i = 2. Assuming statistical independence this probability is time invariant but for 
the cellular automaton the situation is quite different. 

If two cars meet at a site there exist two possibilities: either both cars continue, or one 
of them continues and the other one has to give way. There is no chance that both cars 
meet again one time step later. In the extreme case of only two cars, i.e. N = 2, the 
probability for them to meet by random distribution is thus effectively reduced by a factor 
of two in comparison with the automaton dynamics. Given only a small number N of cars, 
i.e 77 = N/M <C 1, and provided only two of them p and q meet occasionally one has to 
consider in the following time step a system where only N — 2 cars feel the interaction with 
(in our case: can meet occasionally) N — 1 other cars and 2 cars which can meet only N — 2 
cars, i.e. p is not allowed to meet q and q cannot meet p. Hence, the system has some kind 
of memory and therefore it is not ideally Bernoullian. 

Obviously this effect plays a crucial role only for small N since otherwise the relative 
difference between iV — 1 or iV — 2 is negligible. Speaking in terms of density this means 
that the effect is substantial only for very small densities. Then however, clusters larger 
than i = 3 are suppressed by small likelihood. The net effect of these correlations can be 
taken into account through counting all possible configurations which can be realized after 
two cars have met. This procedure is similar to the method employed for the determination 
of the Si (as explained in appendix A). Neglecting contributions related to clusters of size 
larger than three, which is safe for very small densities (see fig 3, explanation below) means 
oii = 1 for i > 4. From simple geometric considerations we find a.\ = 1, a 2 = 1/2 and 
« 3 = 7/4. 

Inserting expression eq.(5) into (1) yields 

- = 1-^ = 1-E^T t k lP (k u ...,k N ) (6) 

i=l fci,...,fcjv 

and performing the sum over all kj (j ^ i) we find 

N J- a- [m 

i=l JV k i= 
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The distribution p(ki) - declared for i = 0,1,..., [N/i\ - can be calculated using the 
inclusion-exclusion principle [50,51]. This derivation is performed in appendix B and the 
result reads 



We see from (7) that the ingredients we actually need are the first moments (Ki) of the cluster 
distribution. They can be derived calculating a generating function for the (descending 
factorial) moments, denoted H.\ M \z,x) [50]. This calculation is performed in appendix C. 
The result reads 



In fig. 3 we have depicted the expectation value of the cluster distribution (Ki) (eq. (9)) 
divided by the system size M = L x L. The solid lines correspond to the function given 
in eq. (9) and the points represent related data taken from our simulation on a 50 x 50 
torus. The values for i = give the probability for a site to be empty. For i = and i = 1 
the analytical and numerical results agree perfectly, for larger i the discrepancies between 
the solid lines and points are a direct consequence of the deviations from the assumed 
independent behavior as explained above. Note that the probabilities in fig. 3 are plotted 
using a logarithmic scale. The observed discrepancies for clusters of size larger than 3 will 
hardly have any influence on the overall behavior of the automaton. For values r\ > rj cr 
the simulation data illustrate a breakdown for the small occupation numbers (except i = 0) 
due to the emergence of jams which are clusters of high order. At the transition point the 
occupation rates abruptly drop by three orders of magnitude. The precise position of this 
transition point is masked by the finite size of N used in the simulation. Increasing the car 
density r\ beyond the blurred critical zone results in a slow rise of the occupation number for 
the small clusters and simultaneously, in a slow decline of the number of empty sites. This 
is due to the increasing number of cars still moving around the jammed crossings. 






(N-i) 



(9) 
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SiOi (N\ 1 / 1 N^-O . . 

E-kr , !-17 • ( 10 ) 



Inserting equation (9) into (7) yields the following formula 

Slcti (N\ 1 A_J_ 
;=1 iV \ J M^" 1 ) V M , 

The last step will be to substitute in (10) the density r\ for the number of cars according to 

N = rjM which results in 

SiOii ( 1 \ 

^ =1 -g^r v'-m)-- 

i - 1\ 1 x^-O 



The main result of our statistical description, namely the mean velocity as a function of the 
car density (eq. (11)), is plotted in fig. 4. The solid line shows the function given by eq. (11) 
and the points are the data taken from the simulation on the 50 x 50 torus (compare fig. 2. 
In the low density region we find a nice agreement. The closer one approaches the transition 
point the more eq. (11) overestimates the simulation results. This can be understood by 
recalling the fact that for densities close to the critical value the effective decomposition of 
the automaton dynamics into the cluster dynamics and the single site dynamics becomes 
inappropriate due to formation and dissolution of short lived jams (like critical fluctuations). 
They tend to increase the occupation number of larger clusters and hence, lead to an effective 
decrease of the average velocity. 

In the thermodynamic limit - which means M, N — > oo and keeping r\ = N/M constant 
- we drop terms of order 0(M~ 1 ) and end up with 

v w 1 - Y] -V V { ' exp(-r/) (12) 
, i /! 

where the exponential corresponds to the dominant contribution of the last factor in (11). 
This formula can be expressed by a series expansion in powers of r\ 



rjM oo / -i \k £ 

v M ~°° i _ y: E . ,, v (t - 1+k) (13) 



i=l fc=0 

( 6202 \ ( $2^2 _ S s a s \ 2 

_v 2 r v_2 e r 

6 2 a 2 S 3 a 3 6 4 aA 3 

-1 — r + ^r) v + ~-- 
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Note that there is no constant term in the double sum since Si = 0. Equation (13) can be 
truncated at a given order of r\ hence, giving rise to an approximation scheme. In fig. 5 we 
plotted the linear (dashed), the quadratic (dotted) and cubic (long dashed) approximations 
together with the full formula (fat solid). The range of reliability visibly increases with 
increasing the order of approximation. On the other hand this scheme clearly illustrates 
that the functional relation between mean velocity and global density definitely is nonlinear^ 
except for very small car densities. 

V. CONCLUSION 

We investigated the behavior of a two-dimensional cellular automaton with periodic 
boundary conditions to simulate traffic flow in cities. The automaton mimics realistic traffic 
rules which apply to our everyday experience of vehicular traffic. In agreement with similar 
models we found a slow decay for the mean velocity (v) as a function of the global traffic 
density r\. At a critical threshold value r\ cr the mean velocity collapses abruptly and the 
system transits into another regime of global behavior which we call the jamming regime. 
Beyond the critical density r\ cr the average velocity is very small and it declines further 
when increasing the density. By simulating different sizes of automata we could exclude the 
influence of finite size effect provided the density value is not rather close to the critical 
density. 

Applying combinatorics and statistical methods for the description of the system in the 
low density regime the analytical calculation performed in section IV yielded the average car 
velocity as a function of the car density - (v) (77) - which nicely agrees with the values of the 
numerical simulations in section III. Since the derived functional relationship between mean 
velocity and global traffic density was based on very general assumptions (conservation of 
the number of cars, weak spatial correlations in the low density regime, ergodicity) we expect 
this description to be valid for a wider class of traffic systems. Note that the specific structure 
of traffic rules enter the description only through the average delays Si. The only nonlocal 
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ingredient arises from the restriction that cars have to stop in case the next desired queue at 
a neighboring site is totally filled up. But this nonlocal character only becomes substantial 
for densities close to the critical value where it causes long range spatial correlations. 

Several authors (e.g. [52]) assert that the average velocity v in the low density regime is a 
linear function of the density 77 and indeed the simulation results seem to support this obser- 
vation. A more detailed analysis however reveals that this function definitely is nonlinear. 
An analytic expansion shows (eq. (13), fig. 5) that the results become dramatically wrong 
when truncating the formula after the linear term. Surprisingly the analytic description 
gives good results even for densities not too small where the nonlocal effects of the dynamics 
and hence, long range correlations spoil the basic assumptions of our treatment. 

For still higher values of the density r\ the results of our simulation agree with results 
reported in literature (e.g. [52] and many others), but no longer with our approximate theory. 
Above a critical density we observe an abrupt transition of the system into the jammed 
regime where the averaged velocity is close to zero. This regime however is beyond the 
scope of our statistical description and has to be investigated starting from other approaches 
(e.g. nucleation processes). 
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APPENDIX A: CALCULATION OF THE AVERAGE NUMBER OF WAITING 
CARS ^7 FOR CLUSTERS OF SIZE / = 1 AND 1 = 2 

In table II we present the average number Si of cars which are prevented from moving 
given i cars meeting at a crossing. Obviously the number Si is confined to the interval [0, i\. 
The average is performed with respect to all possible configurations which are assumed to 
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possess equal probability. The case i = 1 is rather simple: If there is only one car at the 
crossing it has never to give way hence, 8\ = 0. The case i = 2 is not that trivial and the 
computation of 62 requires some combinatorial effort. Since each of the cars can wait at one 
of the four sides of the crossing and can proceed in one of the three directions (left, straight 
on, right) we have 4x3x4x3 = 144 different situations of equal probability. Each of 
these situation can be assigned to one of the situations in the left column of table I. We 
now explain the columns of table I: 

1. A sketch of the configuration given by a related graph. 

2. The number of the different realizations that relate to the graph in the first column 
assuming first that the cars are distinguishable. The first factor A in the form AxBxC 
originates from the symmetry according to rotation by ^7r, 7r and |7r. The second factor 
B denotes the number of choices for the first car (which is now assumed to come from 
South) and the last factor C gives the number of situations possible for the second 
car. 

3. When we take into account that the cars in fact cannot be distinguished the number 
of different events that belong to the figure in the first column has to be multiplied by 
the factor given in the third column. 

4. The total number of events for the whole class denoted by the graph, i.e. (column 2x 
column 3). 

5. The number of cars stopped by the traffic rules at this crossing in one of the possible 
realizations (either one or none). 

6. The total number of cars which have to stop related to the whole class denoted by the 
graph, i.e. (column 4x column 5). 

From the last line in table I one sees that there are 144 different situations of equal 
probability which amounts to 288 cars. Furthermore we see that 92 of those 288 cars have 
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to stop. Hence, the average number of the cars stopped for 2-clusters is S 2 = 92/144 as 0.639. 

For the clusters of size % = 3 the average number of cars which are prevented from moving 
<5 3 can be calculated in analogy however the corresponding table contains about 8 times as 
much situations as table I. Therefore we do not want to present it here. The result is 
<5 3 = 1.479. For cluster of size larger than 3 the corresponding <5j have been calculated by the 
method of complete enumeration using a computer. These empirical results are collected in 
table II. For i = 1 . . . 3 they are identical with the analytical results. 

For completeness sake we remark that for large cluster sizes (i > 8) there exists a 
negligible probability that one of the directions of the crossing is not occupied. Typically 
there are many cars waiting in each of the four queues. In those very likely cases our traffic 
rules allow exactly one (randomly chosen) cars to drive, all other i — 1 cars have to stop. 
Therefore we find Si i — 1. Certainly the contributions of such big clusters will only play 
a minor role in our calculations due to very small likelihood. 

APPENDIX B: THE DERIVATION OF THE CLUSTER PROBABILITIES 

The probability p (ki) to find exactly k sites which are occupied each by i (independently 
moving) cars in a system of M crossings and N cars is given in eq. (8). The derivation of this 
formula employing the inclusion-exclusion principle will be performed in this appendix. Note 
that the problem to find p (ki) is different from the trivial problem to find the probability 
for the ki crossings which are occupied by at least i cars. 

The inclusion-exclusion principle relates the probabilities for a finite number of sections 
of events to the probabilities of an exact number of events. 

Let (Q, P) be a probability space, A\, . . . , A N be events, and for arbitrary {mi, . . . , nij} C 
{1, . . . , iV} let P(A mi fl . . . fl A m .) be known probabilities. We define 

Sj= P(A mi n...nA mj ) (Bl) 
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Moreover, let be B n = {u e Q : co e A m for exactly n values of A:}. Then the inclusion- 
exclusion principle asserts that 

P(B n ) = £ (-l)W-") (B2) 

Applied to our case the Sj are the probabilities that j cells each contain i cars and the 
rest of the cars, i.e. (A — ji) cars, are distributed arbitrarily among the remaining (M — j) 
cells; this probability can be derived with ease and reads 

Sj = (/!)•' I.Y .,/)! (m) ( X " m) (B3) 
To explain this formula we first index each of the A cars with the numbers 1, . . . , A which 
yields the factor AH. Since the i cars within each of the j cells are not distinguishable we have 
to divide this factor by (z!)- 7 and since the remaining (A — ji) cars are not distinguishable 
too additionally by (A — ji)\. After having numbered the cars we subsequently fill the first 
cell with cars 1, . . . , z, the second cell with cars (i + 1), . . . , (2i) and so on. In this way we 
distribute the cars numbers 1, . . . , ji among the urns 1, . . . , j yielding the factor M'- - - 7 *-'. The 
remaining (A — ji) cars are distributed successively among the remaining M — j cells at 
random which explains the factor [(M — j) / M]^^^. 

Insertion of this probability into the inclusion-exclusion formula (B2) yields 

M / i \ 

j — ki 

M N W W (N-ji)\ [ ] 

Because of the generalized definition of factorials - using the gamma function - the sum 
effectively only ranges from ki up to |_^AJ which is sensible. 

APPENDIX C: CALCULATION OF THE FIRST MOMENTS OF THE CLUSTER 

DISTRIBUTION {K^ 



In the following we will derive the mean value of the cluster distribution (Ki) which is 
given in eq. (9). 

19 



The generating function for the (descending) factorial moments of the cluster distribution 



we use is 



oo oo 



N 



i*rw) = £ e 

N=0 ki=0 



M N z 



x hi p(ki,N) 



(CI) 



where p(ki, N) is the probability to find a value hi for the stochastic variable Ki when trying 
with iV balls. The benefit of such a complicated looking generating function is that the sums 
can be performed yielding an analytical expression namely 

M 



exp z +— (x - 1) 



(C2) 



For the explicit derivation the reader is referred to the book by Johnson and Kotz 3 
p. 116ff [50] . Moreover, this function is related to the (descending) factorial moments 
of the stochastic variable Ki according to 



AH d r 



(H< M \x,z) 



x=l 



M N dx r 

OO oo 

= E zN E ^ -!)••• ( k i ~ r + !) p(^ N ) 

N=0 ki=0 
oo 

= E^ (Hki-l)... (ki-r + 1)) 
hence, the first moment (Ki) is given as the coefficient of z N in 



(C3) 



N\ d f „ {M)i 



M N dx 



i-(Hr\z,x) 



x=l 



(C4) 



Consequently we calculate 



N\ d 



M N dx 



exp z +—(x - 1) 



M 



x=l 



exp 



[z(M-l) 



^ i\ k\ M^- 1 ) 



(C5) 



3 one has to identify j = i,m = M,n, = N,Mi = K t ,g = ki,Pi[Mj = g\n] = p(ki) 
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and finally arrive at 

AH (M - l)^" 1 ) 



o: \iy — l ): ivi >■-• 
/N\ 1 / 1 \ ( N -^ 

(C6) 



- 



21 



REFERENCES 

[1] J. Treiterer. Techn. Report PB246, Columbus, Ohio, 1975. 
[2] T. Musha and H. Higuchi. Jap. J. Appl. Phys., 15:1271, 1976. 

[3] A. D. May. Traffic Flow Fundamentals. Prentice Hall, Eaglewood Cliffs, N. Y., 1993. 

[4] W. Leutzbach. Introduction to the theory of traffic flow. Springer Verlag, Berlin, 1988. 

[5] M. J. Lighthill and G. B. Whitham. Proc. Roy. Soc. A, 229:281, 1955. 

[6] M. J. Lighthill and G. B. Whitham. Proc. Roy. Soc. A, 229:317, 1955. 

[7] I. Prigogine and R. Herman. Kinetic Theory of Vehicular Traffic. Elsevier, New York, 
1971. 

[8] R. Herman and I. Prigogine. Science, 204:148, 1979. 

[9] S. Wolfram. Theory and Applications of Cellular Automata. Elsevier, New York, 1971. 
[10] K. Nagel and M. Schreckenberg. J. Physique I, 2:2221, 1992. 

[11] R. Wiedemann. Simulation des Strafienverkehrsflusses, volume 8. Schriftenreihe des Ins- 
tituts fur Verkehrswesen, Universitat Karlsruhe, 1974. 

[12] T. Nagatani. J. Phys. Soc. Jap., 62:3837, 1993. 

[13] K. Nagel and H. J. Herrmann. Physica A, 199:254, 1993. 

[14] K. Nagel. J. Mod. Phys. C, 5:567, 1994. 

[15] K. Nagel. Transportation Research Board, 94 09 01, 1994. 

[16] K. Nagel and A. Schleicher. Appl. Math, and Comp. Sc. (in press), 1993. 

[17] K. Nagel and S. Rasmussen. In R. A. R. A. Brooks and P. Maes, editors, Artificial Life 
IV: Proceedings of the Fourth International Workshop on the Synthesis and Simulation 
of Living Systems, Cambridge, 1994, P a g e 222. MIT Press, Cambridge, MA, 1994. 

22 



[18] A. Schadschneider and M. Schreckenberg. J. Phys. A, 26:L679, 1993. 
[19] Z. Csahok and T. Vicsek. J. Phys. A, 27:L591, 1994. 
[20] J. Walker. Scient. Am., 247:206, 1982. 

[21] B. Sick. Dynamische Effekte bei nichtlinearen Wellengleichungen mit Anwendung in der 
Verkehrsflufitheorie. PhD thesis, Universitat Ulm, 1988. 

[22] B. S. Kerner and P. Konhauser. Phys. Rev. E., 48:2335, 1993. 

[23] E. Ben-Nairn, P. L. Krapivsky, and S. Redner. PRE, 50:822, 1994. 

[24] D. Helbing. Complex Systems, 6:391, 1992. 

[25] D. Helbing. preprint, 1994. 

[26] R. Kiihne. Physik in unserer Zeit, 15:84, 1984. 

[27] R. Kiihne. Gunnoszillator, Fernstrafienverkehr und Flachwasserwellen. In 
Fruhjahrstagung der Deutschen Physikalischen Gesellschaft, Fachausschuss Dynamik 
und Statistische Mechanik, 1988. 

[28] M. Takayasu and H. Takayasu. Fractals, 1(4):860, 1993. 

[29] Y. Sheffi. Urban Transportation Networks: Equilibrium Analysis with Mathematical Pro- 
gramming Methods. Prentice-Hall, New Jersey, 1985. 

[30] H. S. Mahmassani, R. Jayakrishnan, and R. Herman. Transpn. Res. A, 24:149, 1990. 

[31] M. Hilliges, R. Reiner, and W. Weidlich. In: Alain Pave, editor, Modelling and simula- 
tion ESM 93, p. 505, Lyon, 1993. European Simulation Multiconference. 

[32] R. Haberman. Mathematical Models: Mechanical Vibrations, Population Dynamics, and 
Traffic Flow. Prentice-Hall, New Jersey, 1977. 

[33] T. Nagatani. J. Phys. A, 26:6625, 1993. 



23 



[34] M. Bando, K. Hasebe, A. Nakayama, A. Shibata, and Y. Sugiyama. Phys. Rev. E., 
51:1035, 1995. 

[35] O. Biham, A. A. Middleton, and D. A. Levine. Phys. Rev. A, 46:R6124, 1992. 

[36] S. Tadaki and M. Kikuchi. Phys. Rev. E, 50:4564, 1994. 

[37] K. H. Chung, P. M. Hui, and G. Q. Gu. Phys. Rev. E., 51:772, 1995. 

[38] T. Nagatani. J. Phys. Soc. Jap., 62:1085, 1993. 

[39] M. Fukui and Y. Ishibashi. J. Phys. Soc. Jap., 62:3841, 1993. 

[40] H. F Chau, P. M. Hui, and Y. F. Woo. preprint IASSNS-HEP-95/07, 1995. 

[41] T. Nagatani and T. Seno. Physica A, 207:574, 1994. 

[42] J. A. Cuesta, F. C. Martinez, J. Molera, and A. Sanchez. Phys. Rev. E., 48:R4175, 1993. 
[43] J. M. Molera, F. C. Martinez, J. Cuesta, and R. Brito. Phys. Rev. E., 51:175, 1995. 
[44] G. Peng and H. J. Herrmann. Phys. Rev. E., 49:1796, 1994. 
[45] G. Peng and H. J. Herrmann. Phys. Rev. E., 51:1745, 1995. 
[46] T. Poschel. J. Physique, 4:499, 1994. 

[47] K. Nagel. High-speed Micro simulations of Traffic Flow,. PhD thesis, Universitat Koln, 
1995. 

[48] D. C. Hong, S. Yue, J. K. Rudra, M. Y. Choi, and Y. W. Kim. Phys. Rev. E., 50:4123, 
1994. 

[49] R. v. Mises. Rev. Fac. Sci. Univ. Istanbul, 4:145, 1939. reprinted in Selected Papers of 
Richard von Mises, Vol. 2. Providence, R.I.: American Mathematical Society, 313 1964. 

[50] J. N. Johnson and S. Kotz. Urn Models and Their Application. Wiley & Sons, New 
York, London, Sydney, Toronto, 1977. 

24 



[51] U. Krengel. Einfuhrung in die Wahrscheinlichkeitstheorie und Statistik. Vieweg Verlag, 
Braunschweig, second edition, 1990. 

[52] F. C. Martinez, J. Cuesta, J. M. Molera, and R. Brito. Phys. Rev. E., 51:R835, 1995. 



25 



TABLES 

TABLE I. All situations which might occur when two cars meet at a crossing, their frequency, 
their symmetry according to rotation and due to indistinguishability of the cars and the number 



of c 


ai 


"s which have to stop in the current situation (explanation in the 


text). 




graph rot. symm. 


dist. ^events 


waiting 


^stopping 


















4x3x3 


1 36 


1 


36 


- — 




4x1x2 


2 16 


1 


16 






4x1x1 


2 8 












4x1x3 

A 


2 24 


1 


24 


V 


^ 4x1x3 


2 24 










^ 2x1x1 


2 4 












4x1x2 


2 16 


1 


16 






2x1x1 


2 4 












4x1x1 


2 8 










2x1x1 


2 4 










total 


144 




92 
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TABLE II. The average number of cars which are stopped Si when i cars meet at a crossing. 



^cars meeting i Si 



1 


0.0 


2 


0.638889 


3 


1.479167 


4 


2.467014 


5 


3.524740 


6 


4.604709 


7 


5.683838 


8 


6.752964 
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FIGURES 
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FIG. 1. The schematic plot of a crossing 
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FIG. 2. The simulation result: The mean velocity v sim of iV cars on a periodic lattice of size 
IxLvs the global traffic density rj = N/ {L x L). Except for values around the critical density rj cr 
the curve is rather smooth. In the vicinity of the transition zone critical fluctuations, i.e. temporary 
jams, cause an irregular relationship. This can be understood from the fact that the typical lifetime 
of jams becomes of the same order as the simulation time. 
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FIG. 3. The average number of clusters of size i normalized to the lattice size M = L x L 
as a function of the global traffic density rj = N/M. The lines represent eq. (9) while the points 
correspond to numerical data achieved by a simulation on a 50 x 50 lattice (torus). The curve i = 
gives the probability for a site to be empty. 
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0.0 1.0 2.0 

FIG. 4. The average velocity vs the global traffic density for the low density regime (non 
jamming traffic). The curve displays the result from equation (11), the points show data taken 
from the simulation. 
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FIG. 5. The average velocity as derived from the statistical description (eq. (9)) (fat line) and 
truncated series expansions: linear (dashed), quadratic (dotted) and cubic (long dashed). This 
plot clearly illustrates a nonlinear relationship between average velocity and global density. 
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